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Abstract 

We present a theoretical study on the relation between the size of the rare earth ions, often 
known as chemical pressure, and the stability of the coherent Jahn- Teller distortions in undoped 
perovskite manganites. Using a Keating model expressed in terms of atomic scale symmetry modes, 
we show that there exists a coupling between the uniform shear distortion and the staggered 
buckling distortion within the Jahn- Teller energy term. It is found that this coupling provides a 
mechanism by which the coherent Jahn- Teller distortion is more stabilized by smaller rare earth 
ions. We analyze the appearance of the uniform shear distortion below the Jahn- Teller ordering 
temperature, estimate the Jahn- Teller ordering temperature and its variation between NdMnO,3 
and LaMn03, and obtain the relations between distortions. We find good agreement between 
theoretical results and experimental data. 

PACS numbers: 75.47.Gk, 74.62.Dh, 64.70.K-, 61.50.Ks 
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I. INTRODUCTION 



Since the discovery of the colossal magnetoresistance effect, a lot of attention has focused 
on a class of materials known as perovskite manganites.- - - These materials have the chemical 
formula in the form of REi^ x AK x Mn03, where RE and AK represent the rare earth and 
alkali metal elements, and have a perovskite structure. One of the major research themes for 
these materials is the relation between their physical properties and the average size of ions 
at the RE / AK site, often known as the chemical pressure effect. The size of the RE / AK 
ion is usually parameterized by a tolerance factor and one of the most important phase 
diagrams for these materials has been the one in the temperature versus tolerance factor 
plane for a fixed 30% (x = 0.3) doping ratio.- The RE/ AK ions with size smaller than the 
space created by the surrounding MnC>6 octahedra induce buckling of the Mn-O-Mn bonds, 
observed through various structural refinement analyses. 

How this buckling distortion affects the properties of manganites has been controversial. 
It is well known from experimental observations that there is strong competition between 
the insulating phase with a coherent Jahn- Teller (JT) distortion and the metallic phase 
without such distortion.- So far, most attention has centered on the impact of the buckling 
on the metallic phase, in particular, the possible change in the effective Mn-O-Mn electron 
hopping parameter and the band width. 7 However, there has been a debate on whether 
the change of the hopping parameter due to the Mn-O-Mn bond angle change of several 
degrees would be significant enough to explain the observed metal-insulator transition.-"— 
A less studied effect of the Mn-O-Mn bond buckling, except for a few early efforts based on 
experimental data,— is the possibility that the buckling distortion may significantly stabilize 
the insulating phase with a coherent JT distortion. The main goal of this paper is to examine 
such a possibility with a simplified model of the perovskite manganites. To be specific, we 
analyze the interplay between the JT ordering and chemical pressure for undoped perovskite 
manganites. The study on undoped manganites is merited because they are not only parent 
compounds of doped perovskite manganites, but also because one of the first multiferroic 
materials discovered is an undoped manganite, TbMnOs, with a relatively small rare earth 
element.— Therefore, the chemical pressure effect in undoped manganites reported in this 
paper would also be relevant for future studies on how the multiferroic property appears in 
REMn0 3 with small RE ions, as well as for the effect of chemical pressure on the distorted 
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FIG. 1. Two-dimensional model for the perovskite structure considered in the text, 
insulating phase of doped manganites. 

II. MODEL SYSTEM AND ENERGY EXPRESSION 

We study a two-dimensional (2D) model for the perovskite structure which incorporates 
both buckling and the JT distortions. We define a 2D perovskite structure shown in Fig. CD 
which includes the following aspects of the 3D perovskite structure for undoped manganites: 
(1) symmetry breaking distortion of oxygen ions around Mn ion, (2) chemical pressure effect, 
which is the attraction of surrounding oxygen ions toward the small rare earth ions, and (3) 
the rotation of oxygen ions with alternating directions around Mn ions, which is, in effect, 
the buckling of Mn-O-Mn bonds. 

We apply the recently developed atomic scale description of lattice distortions^*^ to 
describe the elastic energy of the system. In this approach, atomic scale modes of lattice 
distortions and their constraints are used instead of displacement variables. The structural 
motifs can be chosen in any convenient way as long as they have the symmetry of the 
crystal structure. We choose two "structural motifs" , shown in Fig. |2) one consists of one 
Mn ion and four surrounding O ions and the other comprises of one rare earth (RE) ion 
and four surrounding O ions. We obtain seven symmetry modes for each motif, shown in 
Fig. [3] for the MnC>4 motif. 1 - Similar symmetry modes are defined for the REO^ motif and 
are distinguished with primes on the symbols in this paper. The modes defined for each 
plaquette on the lattice are constrained by each other because neighboring motifs share 
ions, which leads to constraint equations between the Fourier components of the modes. In 
terms of these fourteen modes and constraint equations, any distortion of the 2D perovskite 

3 



(a) O (b) O 

O • O O O ° 

O O 
around Mn around RE 

FIG. 2. Two structural motifs chosen for the 2D structure shown in Fig. [TJ 
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FIG. 3. Distortion modes for the motif around the Mn ion in Fig. [2j Similar distortion modes, e'i, 
e' 2 , e' 3 , s' x , s' y , t' x , t' y are defined for the motif around the RE ion. 



structure shown in Fig. [T] can be described. 

For the current study, since we are interested in the ordered state, we consider distortions 
with wavevectors k = (0, 0) and (n, it) only. For these wavevectors, the constraint equations 
are as follows, where we use subscripts and s to represent k = (0,0) and (tt, tt), respectively: 

e 10 = e 10 ' e 20 = e 20> e 30 = e 30> S x0 = ~ S xOi S y0 = ~ S yCn e ls = ~ e 2s5 e ls = — e 2s; s xs = s ys = 

s' xs = s' ys = 0. Rest of the modes are unconstrained, particularly, e 3s and e' 3s . We search 
for the interplay between the staggered deviatoric distortion mode e 3s and the staggered 
rotation of O ions around Mn ion (or equivalently staggered Mn-O-Mn bond buckling mode) 
e' 3s , where the latter is due to the compression e' 10 = eio by small rare earth ions. Therefore, 
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FIG. 4. Four distortion modes considered in the current study: (a) uniform dilatation mode eio, 
(b) uniform shear mode e 2 o, (c) staggered deviatoric mode e^ s , (d) staggered buckling mode e' 3s . 
All figures are drawn for the positive values of the modes with the Mn site at the left bottom corner 
chosen as the origin. 



we limit ourselves to the modes e' 10 = e w , e 20 = e 20 , e 3s , and e' 3s , shown in Fig. HJ We 
include the uniform shear mode e 2 o = e 20 because it is coupled to e' 3s through the JT term, 
as will be discussed later in this paper. 

Even though it is possible to analyze an energy expression including higher order 
symmetry-allowed anharmonic energy terms, such a method would generate many parame- 
ters and would make the model less predictive. Therefore, we start with a Keating model 
with a small number of parameters,— 1 ^ and map the Keating model onto the approach 
based on the symmetry modes. In the Keating approach, the elastic energy is represented 
in terms of bond length and bond angle changes from equilibrium. For our 2D perovskite 
structure, we consider the following set of Keating variables and harmonic moduli for each 
Mn ion, as shown in Fig. 5l n (n — 1, 2, 3, 4) and modulus a\ for Mn-0 bond length 
change, S9 n (n = 1, 2, 3, 4) and &i/4 for 90° O-Mn-0 bond angle change, 5r n (n = 1, 2, 3, 
4) and a 2 for RE-0 bond length change, and 5ip n (n = 1, 2) and 6 2 /4 for 180° Mn-O-Mn 
bond angle change. We note that the MnC>4 motif is considered as relatively stiff compared 
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FIG. 5. The Keating variables considered for each Mn ion. l±, I2, h, and I4 represent the Mn-0 
bond lengths. 9±, 62, #3, and #4 represent the O-Mn-0 bond angles. r±, r2, r%, and represent 
RE-0 bond lengths. ip\ and if2 indicate Mn-O-Mn bond angles. 

to other components of the structure, so that a\ ^> a 2 and bi ^> 6 2 . 

We consider the following Keating elastic energy expression per Mn ion, 

^elastic =^1 E (^n) 2 + ^l E ( 66 n/2) 2 
^ n=l, 2,3,4 ^ n=l, 2,3,4 



+ iU E (^n) 2 + ^2 E i^n/2) 2 . 
/ n=l, 2,3,4 / n=l,2 



(1) 



We express the Keating variables in terms of eio, e2o, ^ s and e' 3s . For example, we obtain 

(2) 



5h = 



50i = tan" 1 



1 + eio + e 3s ) 2 + (e 20 + e 3s ) 2 - 1 



e 2 o + e' 



3s 



+ tan 1 



e 20 



'3s 



(3) 



l + eio + e 3s y V 1 + e io _e 3s / 

The Taylor expansion of -Elastic m terms of e 10 , e 2 o, e 3s and e' 3s produces all the terms of 
any order. We make an approximation that 62 is much smaller than other parameters, as 
mentioned above, and drop the terms with b 2 . We keep all harmonic order terms and select 
the cubic and quartic order terms that are responsible for the Mn-O-Mn bond buckling 
instability, which are shown below as E^, E cu u c , and Quartic- 

We further define the JT energy per Mn ion E JT and the energy associated with the 
tolerance factor per Mn ion E to \ as follows: 

A, 



Ejt — 



\5l\ + 5/3 — 5I2 — SI4 1, 



(4) 



E to \ = -jiSn + 5r 2 + Sr 3 + <5r 4 ), 



where we define "chemical pressure" 



p = C[(l-t). 



(5) 



(6) 



The parameter t is a two-dimensional analog of the tolerance factor for the 3D perovskite 
structure, and the coefficient C[ represents the coupling between the average RE-0 bond 
length and the tolerance factor t. The "chemical pressure" p induces the shortening of the 
average RE-0 bond length due to small RE ions. We also define the JT distortion mode 



e JT = (5h + 5l 3 - Sl 2 - SU)/2, 



(7) 



which represents the anisotropic bond length change, similar to the 3D JT distortion modes, 
often written as Q 2 and Q 3 (Ref. |201 ). The expression E^ = — A|ejT | is based on the 



3D JT energy -Ejt,3d = — Ql + Ql which is obtained after minimizing JT electron- 
lattice coupling energy in manganites with respect to_the e g orbital state.— In undoped 



manganites, Q3/Q2 is about 0.3 - 0.4 (Refs. 



21 



and 



22), 



which allows an approximation 



£jt,3D ~ -Aq|Q 2 |[1 + (<5s/<52) 2 / 2 ]- Further neglecting the small (Q 3 /<52) 2 / 2 term, we get 
the 2D analog of the JT energy £jt, with the 2D JT distortion ejx corresponding to 3D JT 
distortion Q 2 except for a normalization factor difference. 

We expand Ej T and E to i in the form of a Taylor series in e w , e 2 o, e 3s and e' 3s . Only the 
leading order energy terms being kept, our total energy expression per Mn ion E tot is given 
below. 



Etot 



-^har + -EjT + Ei 



tol 



E 



cubic 



+ E 



quartic , 



1 1 

-(ai + a 2 )(eio) 2 + -(46 1 )(e 20 ) 2 



-Ejt = -A|e 3s + e 2 oeg s 
E to \ = peio, 
-Ecubk = 2 aiei o( e 3s) 2 ' 



E 



quartic 



1 ax, , 



4 2 



"3s ) 1 



(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



where the relation 



e JT ~ e 3s + e 20 e' 3s 



(14) 



is used for Ejt. The physical origin of the coupling between e2o and e' 3s is important for the 
current study and is explained in more detail in Sect. IV Al 



III. ESTIMATION OF PARAMETERS 



In this section, we present our estimate of the parameters. We choose the Mn-Mn distance 
before the distortion, which is around u — 4 A, as 1. Therefore, ei , e 2 o, e 3s , and e' 3s are 
unitless, and a%, bi, ci2, &2, and A have the unit of energy. The parameter a\ can be estimated 
from the Mn-0 bond stretching phonon mode energy, which is about 70 meV from optical 



measurements.— From fry 2a 1 /mo = 70 meV with mo the mass of the oxygen ion, we obtain 
a\ « 150 eV. We estimate b\ from the elastic modulus, C44. From Ref. 23], C44 ~ 55 - 60 GPa. 
The uniform shear mode e2o corresponds to the conventional e xy /2 (Ref. \24 ). Using the 
identity 1 GPa A 3 = 6.3 meV, we find bi ~ 20 - 25 eV. To estimate 62, we use the results^ 
for ReO"3, which has no RE/AK ion and, therefore, a 2 = and the buckling of Re-O-Re bond 



depends only on 62- According to the analysis in Ref. |25j, the oxygen oscillation along Re-O- 
Re direction has the angular frequency = 905 cm -1 , whereas the oscillation perpendicular 
to Re-O-Re direction has the angular frequency u y Q = 30 cm -1 , from which we can estimate 
62/ai = { UJ o/ UJ 'o) 2 /'^ ~ 0-5 x 10~ 3 . We can expect a similar order of magnitude for 62 in 
manganites, order of 10 _3 ai, for example, 0.2 eV, which is negligible compared to other 
parameter values and justifies neglecting the terms with 6 2 as mentioned above. Various 
probes, such as neutron or optical spectroscopy, indicate the buckling mode frequency in 
manganites of about 35 - 50 meV (Ref. 



261 ). From the analysis of (n, 71) phonon mode for 



our model, we obtain the frequency of buckling mode Wbk = y (2ct2 + 46 2 )/'Tio- Therefore, 
we obtain a<i ~ 30 - 80 eV. For the estimation of A, we match the JT energy gain for our 2D 
model with that for the 3D model to ensure that our 2D model represents the energy scale 
of the 3D materials correctly. For our 2D model AEj^ = — A 2 /(2ai). For the 3D model in 



Ref. 



20i AE JT w —0.29 eV, and therefore, we obtain A w 10 eV. 



IV. INTERPLAY BETWEEN Mn - O - Mn BOND BUCKLING AND THE JAHN- 
TELLER DISTORTIONS 



A. Buckling instability without the Jahn- Teller term 

We find the condition for the buckling instability without the effect of the JT energy 
term Ej T . We take a perturbative approach rather than try to solve high order polynomial 
equations. By minimizing + E to \, we obtain 

(e 10 ) mm '* = V — , (15) 

a\ + a 2 

where the superscript * indicates that the JT term is not yet taken into consideration. 
This isotropic compression of the Mn0 4 motif renormalizes the coefficient of the (e^) 2 
term through the -E cu bic term. From this, we obtain the critical condition for the buckling 
instability, 

p* c = —(a 1 + a 2 ), (16) 
(eio)r'* = --■ (17) 

If p > p*, Mn-O-Mn bond buckling occurs and the quartic order term, -E qua rtic, should be 
considered for the equilibrium e' 3s , 



l(4) min '*l = J— -y/F^c (18) 

V (L\ + 0,2 



V2 V (e 10 )f^ - (e 10 ) min <*. (19) 



The minimized E tot without the Ejt term is given by 

2{ai + a 2 ) 2 \ai + a 2 ai/ 
B. Buckling instability with the Jahn- Teller term 

We now examine how the JT energy term alters the buckling instability. From 
-Ehar + E toX + £ JT , we obtain 

(eio) min = (21) 

ai + a 2 

(e 3s ) min = -, (22) 

CLi 

9 



where we consider the (e3 S ) min > case only. The buckling instability is found from the 
second order terms in e 2 o and e' 3a in E tot : 

^(46i)(e 20 ) 2 + \[a 2 + a 1 (e 10 ) min ](eL) 2 - Ae 2Q e 3s , (23) 

where we assumed (e3 S ) min + e 2 oe' 3s > 0. From the condition 4&i[a 2 + ai(eio) mm ] < A 2 , we 
obtain the critical condition 

p c = — (ai + a 2 ) - — (ai + a 2 ) (24) 
ai 4oiai 



and the buckling distortion occurs for p > p c . Comparing with p* in Eq. ( 116JL we find that 
the JT energy makes buckling more likely. After this buckling instability, we should include 
the -Equartic term to find the equilibrium result. For this, we first minimize E tot with respect 
to the shear distortion e 2 o to obtain 

(e 2 o) mm = ^4- (25) 

Inserting this back, we get an energy expression for E tot only in terms of e' 3s , which gives 
the equilibrium buckling distortion and the minimum energy, 



&s) Umi =\-—-^p-Pc, (26) 
y ai -+- a 2 



V 1 X 1 

j^mii\ £_ 

tot ~ 2( fll + a 2 ) 2a, 



-*(-E—* + *) . (27) 
2 y a, + a 2 ai 46iai / 

Therefore, the energy gain due to the JT energy term is given by 

AE ^ (P-P*J X2 (28) 

JT 2ai 4(a 1 + a 2 )6 1 1 ; 

up to order A 2 . The second term corresponds to the part of AEj^ which depends on the 

size of rare earth ion, or p. This result shows that the small rare earth ion, or large chemical 

pressure, stabilizes the JT distortion. 



V. COMPARISON WITH EXPERIMENTS 



We make comparisons between our model and experimental results. In Sect. IV A\ we 
explain the simultaneous appearance of the uniform shear distortion and the long range JT 
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distortion observed in undoped manganites. 21 In Sect. IV B| we estimate the changes in the JT 
ordering temperature Tjt between LaMn03 and NdMnC>3, and compare with experiments. 
In Sect. IV CI we calculate the ratios between different distortion modes and compare with 
experimental data for LaMnOs and NdMnOa. 



A. Appearance of uniform shear distortion below the Jahn- Teller ordering tem- 
perature 



Experimental data in Refs. |21| andl22l show that the difference between the lattice constant 
a and b along the diagonal directions in the plane appears simultaneously with the long range 
JT distortion below Tjt for both LaMnOa and NdMnOa. This distortion corresponds to the 
uniform shear distortion in our model, related by e 20 = (b — a)/(2v / 2w) with u = 4 A. We 
analyze the coupling between the JT distortion and the uniform shear distortion, which is 
important for the stabilization of JT ordered state by the chemical pressure. In our model, 
such coupling originates from the term e2oe' 3s m e JT m Eq. (fT4"l) or in in Eq. ffTOl) , which 
can be understood as follows. We consider applying a positive €20 shear distortion to the 
lattice, as shown in Fig. [6] by the axis of elongation and compression along 45° and 135°, 
respectively. Such uniform shear distortion makes the Mn-0 bond lengths either longer 
or shorter depending on whether the direction of the bond is closer to the orientation of 
elongation (45°) or compression (135°), except for the bonds with directions right between 
the two directions. If the system does not have (it, it) buckling, as shown in the thin solid 
lines in Fig. El all Mn-0 bonds make equal angles from the axis of elogation/compression, 
and therefore e2o shear distortion keeps all Mn-0 bond lengths equal. This implies that e2o 
distortion alone does not contribute to the JT distortion or JT energy gain. In contrast, 
if the system has a buckling distortion e' 3s with a wave vector k = (tt,tt), as shown in the 
thick solid lines in Fig. [6j the e2o shear distortion elongates Mn-0 bonds marked with / and 
shortens Mn-0 bonds marked with s, depending on whether the bond direction is closer 
to the axis of elongation or the axis of compression, which results in the JT distortion ejT 
with a wave vector k = (tt, tt). If this extra JT distortion is in the same [opposite] phase as 
[to] the deviatoric e^s distortion, in other words, if e2oe 3s and e^ s have the same [opposite] 
sign, this extra JT distortion increases [decreases] the net JT distortion, which explains 



the expression of in Eq. (TT0|) or ejT in Eq. ffT4|) . We emphasize here that the extra 
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Axis of compression 

FIG. 6. Superposition of (ir, w) buckling e' 3s and uniform shear distortion e2o effectively generates 
the extra (n, it) JT distortion of Mn-0 bond lengths, as indicated by the s and I for the shortened 
and elongated bonds, which is responsible for the &20^ s coupling within ejT and the Jahn- Teller 
coupling -Ejt- In the (vr,7r) JT ordered state, this adds up to the e% s deviatoric mode if es s and 
e2oe' 3s have the same sign. This mechanism is responsible for the appearance of the uniform shear 
distortion below the JT ordering temperature, as explained in Sect IV Al The extra JT energy gain 
for the buckled lattice is responsible for the increase in Tjt in REM11O3 with small RE ions, as 
explained in Sect IV Bl 



JT energy gain occurs only when the e2, e% and e' 3 distortions are in the right phase with 
respect to each other. Experiments 21 show that the (n, ir) Mn-O-Mn bond buckling persists 
even above Tjt without much change in size. However, above Tjt, the coherent e% distortion 
does not exist, and therefore the extra JT distortion due to the uniform e2o distortion in the 
presence of staggered buckling distortion would increase the JT energy gain in some regions 
and decrease the JT energy gain in other regions, and does not change the net JT energy. 
In other words, the energy gain due to the cooperative effect between e 3 , e' 3 , and e 2 does 
not exist at T > T JT . We therefore expect that the e 2 o mode does not exist above T JT and 
appears simultaneously with the long range JT ordering, consistent with the experimental 
results. 
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B. Jahn- Teller ordering temperature and its variation between LaMnC>3 and 
NdMn0 3 



It is reported^ that Tjt changes from 750 K for LaMn03 to 1100 K for NdMn03 by 
about ATjt = 350 K. We estimate ATjt from our model to understand how such a drastic 
change of the JT ordering temperature can occur by the increase in chemical pressure. 

We rewrite E^ in Eq. (|2"Tj) for p > p c as follows. 

E min = V 2 a l(p-Pc) 2 (2Q) 

tot 2( ai + a 2 ) 2ai 2(a l + a 2 ) 2 ' 1 ' 

where 

Pc=P* c - Sp c , (30) 
A 2 

Sp c =— — (oi + o 2 ). (31) 
4oiai 

Since 5p c , the change in the critical chemical pressure due to the Ej^ term, is small relative 
to p — p*, with 5p c / {ft — p* c ) ~ 0.3 for parameter values in Sect. IIIH we keep the terms up to 
linear in Sp c only and rewrite according to the origin of each term as follows. 

rpmm zT-min , 77mm , rpmm , rpmm (on\ 

1 v 2 

771mm f /qq\ 

"»"«•" 2a 1 + a 2 ' [66> 
^5? = -~. (34) 

^ = "5 (35) 

^bk,JT,sh = — 7 ; \o — Pc) dpc, 

{di + a 2 ) 

A "' -(P-S)- (36) 



46i(ai + a 2 ) 

The first three terms, T™™ p , ^jt 1 ' an d T^ 111 , represent the energy terms purely due to 
compression, JT distortion, and buckling, respectively. The fourth term is the energy due 
to the coherent buckling, JT and shear distortions, indicated by its dependence on p — p*, 
A and bi, which gives extra stability to the JT ordering due to the chemical pressure. 

To estimate Tjt, we consider a high temperature state with random JT distortions, for 
which the energy can be written in a similar way as Eq. ( 1321) except for the absence of the 
fourth term due to the lack of coherence among distortions as explained in Sect. IV Al 

E to t = E comp + E JT + E hk . {67) 
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We expect E 1 ™ ss and ss T^ 11 since the unit cell vex 



ume and buckling angle 



do not change very much as the temperature crosses Tjt (Ref. l2ll ). Therefore, the energy 
difference between JT ordered and JT disordered state is 

ZTwan 7-imin zriran 7-rmin 7-rmin fOQ\ 
-^tot ~~ ^tot ~ ^JT ~~ ^JT ~~ • C 'bk,JT,sh- 

We first verify that our model gives the correct order of magnitude of Tjt itself. An 
order of magnitude estimate for Tjt can be made from the energy difference between two 
different JT ordered states, one the most favored state and the other relatively unfavored 
state. The most favored state is that with the JT distortion of k = (tt, tt) considered so far 
in this paper and has the JT energy of Ey£ n = — A 2 / (2ai). We choose a state with the same 
size of JT distortion e 3 but with a wave vector k = (0, 0) as a relatively unfavored state, 
with energy Ej^ d = — A 2 /[2(ai + a 2 )]. Using the estimated parameter values, a\ = 150 eV, 
a 2 = 30 - 80 eV, A = 10 eV, we obtain E^f - Ef^ n « 600 - 1300 K, which has the same 
order of magnitude as the experimentally observed Tjt in the range of 750 K - 1100 K. 

For the change in Tjt between LaMn0 3 and NdMn0 3 , the only term in Eq. (|38p 
which changes with the RE ion size is —E^ n JT shcar . Therefore, the JT ordering tem- 
perature variation between LaMn0 3 and NdMnC>3 can be related to — E^ ° T sh (NdMn03) + 
-^bk jt sh(LaMn0 3 ) within a factor of the order of one. We express E^ n JTsh in terms of 



(4J mm > 



^b m k i ,JT,sh = -^[(4) min ] 2 - (39) 

According to the experimental data,— >22»21 the Mn-O-Mn bond angle is 155° for LaMnOs 
and 150° for NdMnC>3, which corresponds to (e 3s ) min of about 0.22 and 0.27, respectively. 
These distortions, along with parameter values A = 10 eV and b\ = 20 -25 eV, result in 
-£bk,jT,sh( N dMn0 3 ) + E^ JT sh (LaMn0 3 ) w 12 - 16 meV w 140 - 190 K. From a classical 
Monte Carlo simulation for the double-well potential model in Ref. |l5|, we find that the 
structural ordering temperature is about twice the energy difference between the distorted 
ground state and undistorted high energy state.— Although such a relation would depend 
on the details of the model, if we assume a similar situation in the current model, the JT 
ordering temperature variation can be estimated as twice the energy difference, therefore, 
T JT (NdMn0 3 ) - T JT (LaMn0 3 ) « 2 x [-^ k in JTiSh (NdMn0 3 ) + £^ T>Bh (LaMn0 3 )] = 280 - 
380 K, which agrees well with the experimental change in Tjt, 350 K. 
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This agreement shows that indeed the JT ordered state is more stabilized when the 
buckling increases for smaller rare earth ions for undoped compounds. The relatively large 
increase in the JT ordering temperature, both in theory and experimental data, shows that 
the interplay between the rare earth ion size and the JT distortion is significant, and should 
be taken into account to explain the well-known temperature-tolerance factor phase diagram 
of both undoped and doped perovskite manganites. 



C. Relation between shear, buckling, and deviatoric distortion 

Equations (J2"2j) and (1251) imply that the following quantities remain constant regardless 
of the variation in chemical pressure: 

(e3s) min = A 

ai 

(e 20 ) min = ± 

( e y min 4&j ' 

(e 20 ) min = Oi_ 
(e 3s ) min (ey min 46i ' 

We calculate these quantities from the experimental data for LaMnOs and NdMnC>3, and 
present the results in Table [H in which we also show the relation between the distortion vari- 
ables in our model and experimental parameters and the estimate of Tjt and ATjt obtained 
in Sect. EH The results show that (e 3s ) min , (e 20 ) mi 7(ey min , and (e 20 ) mi 7[(e 3s ) min (ey min ] 
remain constant within 2 %, 10 %, and 7 %, respectively, in spite of 32 % and 19 % changes 
in (e 2 o) min and (e' 3s ) mm . These values also agree well with theoretical estimates obtained 
from the parameters in Sect JIIIL The results underscore the strong coupling between these 
distortions, in particular, the important role played by the uniform shear distortion in con- 
necting the Jahn- Teller and buckling distortions, an aspect neglected in the literature so 
far. 



(40) 
(41) 
(42) 



VI. CONCLUSION 

From the analysis of a Keating energy expression expanded in terms of the atomic-scale 
symmetry-modes, we find that the effect of small rare earth ion size, known as chemical 
pressure effect, is significant in stabilizing the long range Jahn- Teller distortion in undoped 
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TABLE I. Parameters from experimental data and comparison with theoretical estimates. Exper- 
imental data for the lattice constants, bond lengths, and bond angles for LaMnC>3 and NdMn03 

n n 

are from Ref. 21 and 22, respectively, measured at room temperature. 



Parameters 


LaMn0 3 NdMn0 3 


Theoretical estimates 


Lattice constant, a 
Lattice constant, b 

Long Mn-0 bond length within ab plane, I 
Short Mn-0 bond length within ab plane, s 
Mn-O-Mn bond angle within ab plane, ip 


5.54 A 5.414 A 
5.75 A 5.731 A 
2.718 A 2.20 A 
1.907 A 1.90 A 
155.1° 149.8° 




Mn-Mn distance with ei only, (a + b)/(2y/2) 


3.995 A 3.94 A 


Compressed from u ~ 4 A 


e 20 = (b - a)/(2 v / 2«) 

4s = (f - y)/2 

e2o/e' 3s 

ejT = (I - s)/u 
e^s ~ ej T - e 2 oe' 3s 
e2o/(e 3s e' 3s ) 


0.019 0.028 
0.22 0.27 
0.09 0.10 
0.068 0.075 
0.064 0.067 
1.4 1.5 


A/(46i) = 0.10 - 0.13 

A/oi = 0.067 
ai /(46i) = 1.5 - 1.9 


ATjt 


750 K 1100 K 
350 K 


£unif _ £min = gQQ _ 130 Q R 

-2A£§5 TjSh = 280 - 380 K 



perovskite manganites. We obtain good agreement with the experimental data on the Jahn- 
Teller ordering temperature and the substantial increase of the Jahn- Teller ordering tem- 
perature from LaMnC>3 to NdMnCV We propose that similar effects need to be considered 
to understand the phase diagram for doped perovskite manganites. We also explain the 
appearance of the uniform shear distortion below the Jahn- Teller ordering temperature in 
terms of the coupling between coherent shear, buckling, and deviatoric distortions within 
the Jahn- Teller energy. Moreover, we estimate the ratio between these distortions at low 
temperature, and find good agreement with experimental data for LaMnC>3 and NdMn03, 
which confirms the coupling proposed between them in our model. 
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